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An overview of several recent developments in density functional theory for classical inhomogeneous liquids is 
given. We show how Levy's constrained search method can be used to derive the variational principle that un- 
derlies density functional theory. An advantage of the method is that the Helmholtz free energy as a functional 
of a trial one-body density is given as an explicit expression, without reference to an external potential as is the 
case in the standard Mermin-Evans proof by reductio ad absurdum. We show how to generalize the approach 
in order to express the internal energy as a functional of the one-body density distribution and of the local 
entropy distribution. Here the local chemical potential and the bulk temperature play the role of Lagrange mul- 
tipliers in the Euler-Lagrange equations for minimiziation of the functional. As an explicit approximation for the 
free-energy functional for hard sphere mixtures, the diagrammatic structure of Rosenfeld's fundamental mea- 
sure density functional is laid out. Recent extensions, based on the Kierlik-Rosinbergscalar weight functions, to 
binary and ternary non-additive hard sphere mixtures are described. 
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1. Introduction 

The theoretical study of inhomogeneous classical liquids received a boost through the development of 
classical density functional theory (DFT). Evans' 1979 article 1 1] constitutes a central reference to DPT. The 
proof of the variational principle that underlies the theory, including the existence and the uniqueness 
of the free energy functional, can be viewed as the classical analogue of Mermin's earlier (1965) work on 
quantum systems at finite temperatures This forms a generalization of the Hohenberg-Kohn theorem 
for ground state properties of quantum systems. Historic milestones that predate these developments, 
and that can be re-formulated in classical DFT language, are Onsager's 1949 treatment of the isotropic- 
nematic liquid crystal phase transition in systems of long and thin hard rods 13], and van der Waals' 1893 
theory of the microscopic structure of the liquid-gas interface |4]. 

Both the Hohenberg-Kohn proof and the Mermin-Evans proof start from a variational principle for 
the respective many-body function. In the quantum case, this is the Rayleigh-Ritz inequality for the many- 
body groundstate wave function. In the classical case, the theory rests on the Gibbs inequality for the 
equilibrium many-body phase space distribution. The corresponding functionals are the ground state 
energy in the quantum case and the thermodynamic grand potential in the classical case. Both functionals 
depend (trivially) on the position-dependent external one-body potential. Via an intricate sequence of 
arguments f]],(2fl, which has become textbook knowledge |5], the dependence on the external potential 
is played back to the more useful dependence on the (in general) position-dependent one-body density 
distribution. 
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In 1979 Levy showed that quantum DFT can be obtained in a more compact and straightforward way 
via a method that he called the constrained search |6] (see also |7]). Here, the search for the minimum 
is performed in the space of all trial many-body wave functions. The constraint that is fixed during this 
search is that all trial wave functions considered need to generate the given one-body density distribu- 
tion. Levy's derivation is both rigorous and elegant and constitutes a standard reference for electronic 
structure DFT. Among the impressive number of citations of his 1979 article, there are only very few pa- 
pers that draw connections to classical DFT, being an example, although the method permits rather 
straightforward application to the classical case 19J. See also Percus' general concept of "overcomplete" 
density functional (loll . 

The constrained search method offers two significant benefits over the Mermin-Evans proof. One is 
simplicity, avoiding the reductio ad absurdum chain of arguments. The other is that it yields an explicit 
definition of the Helmholtz free energy density as a functional of the (trial) one-body density distribution. 
This formula, as reproduced in <2.2> below, is explicitly independent of the external potential. The gener- 
alization that facilitates this development is the concept of minimization in the space of all phase space 
distributions under the constraint of a given one-body density i2.4) . The issue of representability of the 
one-body density has been addressed in flT]. It turns out that the simplicity of Levy's method allows one 
to construct more general DFTs. An example is the variational framework developed in |12], which rests 
on the internal energy functional (rather than the Helmholtz free-energy functional), which depends on 
the one-body density and on a local (position-dependent) entropy distribution. A dynamical version of 
this theory, based on linear irreversible thermodynamics and phenomenological reasoning, is proposed 
in (12I]. Several common approximations for free energy functional were transformed to internal en- 
ergy functionals. Having reliable approximations for the functional is a prerequisite for applying DFT to 
realistic (three-dimensional) model systems. The task of constructing usable functionals is different from 
the conceptual work outlined so far. In particular, a calculation of the constrained search expressions 
for the (free-energy or internal-energy) functional would amount to an exact solution of the many-body 
problem. Hence, the importance of these expressions is rather of conceptual nature. 

In 1989 Rosenfeld wrote a remarkable letter in which he proposed an approximate free energy func- 
tional for additive hard sphere mixtures |13]. His theory unified several earlier liquid state theories, 
such as the Percus-Yevick integral equation theory for the bulk structure |5], the scaled-particle theory 
for thermodynamics, and Rosenfeld's own concepts, such as the scaled-field particle theory of JIJ], and 
encapsulated these into what he called fundamental measure theory (FMT). Kierlik and Rosinberg in 
1990 1 14] re-wrote the same functional |15] in an alternative way, using only four (not six, as Rosenfeld) 
weight functions to build weighted densities via convolution with the bare one-body density of each hard 
sphere species. The two strands of FMT were pursued both with significant rigour and effort, see the re- 
cent reviews |16-18]. Rosenfeld's more geometric approach was extended to further weight functions by 
Tarazona |19] in his treatment of freezing. A critical discussion of the properties of the relevant convolu- 
tion kernels can be found in |20]. Very recently, Korden I2ii22j demonstrated the relationship of the FMT 
with the exact virial expansion. 

It is vexing that Rosenfeld himself, who was certainly very versed in the diagrammatic techniques of 
liquid state theory (see e.g. | 23]), apparently neither analysed nor formulated his very own FMT within 
such a framework. The non-local structure of FMT, its coupling of space integrals via convolution, and the 
plentiful appearances of the one-body density distribution(s) seem to constitute an ideal playground for 
formulation in diagrammatic language. A comprehensive understanding of the diagrammatic nature of 
FMT could not help just to ascertain and clarify the nature of the approximations that are involved i2l|, 
but also, and more importantly from a pragmatic point of view, could enable one to construct new func- 
tionals for further model systems. Taking the FMT weight functions as bonds in a diagrammatic formu- 
lation, one immediately faces the combinatorial problems associated with their number, which in the 
Rosenfeld-Tarazona formulation are at least seven per hard sphere species (four scalar, two vector, one 
tensor), which makes the book-keeping task a seemingly daunting one. 

It was shown |24, 25] that one can formulate the Kierlik-Rosinberg form of FMT in a diagrammatic 
way. The concept was applied to one-dimensional hard rods, where it gives Percus' exact result jiiiliTlj . 
as well as to five-dimensional hard hjrpersphere mixtures, where it gives a functional that outperforms 
Percus-Yevick theory, and improves on previous FMT attempts jisll . Two crucial properties of the dia- 
grammatic formulation can be identified. One is the relative simplicity of the diagrams that describe 
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the coupling of the various space integrals in the density functional. The topology of the diagrams is of 
star-like or tree-like shape, hence providing a significant reduction as compared to the complexity of the 
exact virial series. The number of arms is equal to the power in density and the bonds are weight func- 
tions rather than Mayer functions, like they are in the exact virial expansion. The book-keeping problem 
of having to deal with a large number of different weight functions is addressed and rendered almost 
trivial by exploiting the tensorial structure that underlies the Kierlik-Rosinberg form of FMT |[29(], where 
the geometric index of the weight functions is a proper tensorial index, and corresponding (isometric and 
metamorphic) transformations can be formulated |29]. Hence, the fuUy scalar Kierlik-Rosinberg formu- 
lation turned out to be indeed simpler to handle, and easier to generalize for full control of the degree 
of non-locality in the functional. These developments facilitated the generalization of FMT for binary 
non-additive hard spheres |30J to ternary mixtures |31]. 

In the present contribution, we describe the basic ideas underlying the above developments, without 
the full detail that is given in the respective original papers, but with further illustrative examples in 
order to provide an introduction to the subject. The paper is organized as follows. In section |2] Levy's 
method is sketched, and both the free-energy and the internal-energy functionals are defined. A basic 
introduction to the diagrammatic formulation of FMT is given in section [S] including a brief overview 
of applications to non-additive hard sphere fluids in bulk and at interfaces. Section |4] gives concluding 
remarks. 



2. Levy's constrained search in classical DFT 

2.1. Variational principle for the grand potential functional 

Consider a classical many-body system with position coordinates ri,r2,...,rjv and momenta 
pi,p2, . . . ,pAr, where N is the number of particles. Denote the classical trace in the grand ensemble by 



(2.1) 



where h is Planck's constant. For a given total interatomic potential t7(ri,...,rAr) between the particles, 
one defines the intrinsic Helmholtz free energy functional by the explicit expression LgJ 



^ [p] - min 



Tr/ 



( N p2 



(2.2) 



where T is temperature, is the Boltzmann constant, m is the particle mass, p,- - |p, |, and the min- 
imization searches aU many-body (phase space) probabihty distributions /(ri,...,rAr,pi,...,pjv;A^) that 
are normalized according to 

Tr/=1, (2.3) 
and that yield the fixed trial one-body density p(r) via 



JV 



p(r)=Tr/^5(r-ri). 



(2.4) 



i=l 



The relationship 12. 4> is indicated as / ^ p in the notation of d2.2> . It can be shown /i*] that the intrinsic 
free energy functional <2.2> satisfies the Euler-Lagrange equation 



5p(r) 



(2.5) 



Po 



where p is the chemical potential, v{r) is an external one-body potential acting on the system, so that the 
total external potential for a given microstate is Xjlj v{ri), and po(r) is the equilibrium one-body density. 



po(r) =Tr/o ^(5(r-ri), 



(2.6) 
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where the equUibrium many-body distribution /o is given by the Boltzmann distribution, 



Here, the normalization constant is the grand partition sum 

Hn-hN 



(2.7) 



H - Trexp 



(2.8) 



N 



and the Hamiltonian is 

A concise version of the proof of <2.5> is laid out in the following. 



N p2 



(2.9) 



2.2. Sketch of the constrained search proof a la Levy 

It can easily be shown 1 1] that the equilibrium grand potential Oq for Hamiltonian Hn is obtained 
from minimizing Mermin's form of the grand potential functional. 



no = mmTr/ [Hn -fiN+kB Tin/) , 



(2.10) 



where / is a trial many-body distribution. Decompose the right hand side into a double minimization 

no ^ mmminTr f [Hn - IJ.N + kBTlnf) , (2.11) 

P f^p 

where the inner minimization is a search under the constraint that the / generates p via i2A) . For suitable 
HamUtonians as <2.9> with external potential v{r), this can be written as 



Do = minminTrf 



( N p2 N 

E 7^ + + E ^^^i^ -^N+ks T\nf 

\i=l i=l 



In the expression above 



Tr/ 



N 



= Jdr[i'(r)-/i]p(r), 



(2.12) 



(2.13) 



because / ^ p. So we may re-write <2.12> as 



Do^niin-i / dr [t^lr) -/z] p(r) +minTr/ 
P I J f^P 



I N p2 



or 



no = mm IJ dr [ v{r) - ^] p(r) + ^ [p] 



(2.14) 



where ^[p] is given by <2.2) and hence in equiUbrium the Euler-Lagrange equation <2.5) follows. Finally, 
note that from | |2.11> the grand potential density functional can be defined as 



n[p] =min[Tr/(Hjv-pAf+fcBrin/)l. 

f^P 



(2.15) 
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2.3. DFT based on the internal-energy functional 

Using two constraints, rather than one, one can define the internal-energy functional as 



Elp, s] - min < f 



E 

1=1 



2m 



[/(ri,...,riv) 



where the constraint for the density is | |2.4> and that for the local entropy distribution s(r) is 

1 ~ 

s(r) = -A;BTr-^5(r-ri)/ln/. 



The equilibrium value for the local entropy is 



1 



JV 



So(r) = -fcBTr— ^5(r-r,0/oln/o, 

^ i=l 



(2.16) 



(2.17) 



(2.18) 



with the equilibrium many-body distribution /o given in <2.7) . and the total entropy obtained as the spa- 
tial integral So = / drso(r). The Euler-Lagrange equations have the form 



SE[p,s] 



Sp{r) 
SE[p,s] 



^p.- v{r), 



Po.so 



dsir) 



= T. 



(2.19) 
(2.20) 



Po.so 



ExpUcit forms of various internal-energy functionals can be found in 1 12], as can be a dynamic prescrip- 
tion, simUar in spirit to dynamical DFT 13211 . for the joint time evolution of s(r, t], p(r, f), where t is time. 
The proof of (O6)-(2|20) is a straightforward application of Levy's method, as laid out in section [2^ for 
the explicit form see I12n . 



3. Fundamental measure theory for hard sphere mixtures 
3.1. Diagrammatic formulation for additive and non-additive mixtures 

We turn to hard sphere mixtures and use the common splitting of the intrinsic Helmholtz free energy 
into ideal gas and excess (over ideal gas) contributions 

Fl{Pi}] = ksTY, [ drpi(r){ln[pi(r)A^] -1} + Fexc[{p,}], (3.1) 

where p,(r) is the one-body density distribution. A,- is the thermal de Broglie wavelength of species i, 
and Fexc[{Pi}] is the Helmholtz excess free energy functional that is due to interparticle interactions. Its 
lowest order (in density) virial expansion is 

fexc[{p,}] - -- / drdr'pKr)/i;(r-r')p^(r'), (3.2) 

where the Mayer function is defined as fijir) - exp[Uij{r)liksT)] - 1, with Uij{r) being the pair in- 
teraction potential between species / and j. The sum in i3.2\ is over all species For hard spheres 
fijir < aij) - -1 and zero otherwise; here cr,y is the hard core interaction distance between species i 
and;. 

Kierlik and Rosinberg introduced a set of four "weight functions" wo{Ri,r), wi{Ri,r), W2{Ri, r), and 
u>-i{Ri,r), where 7?,- = (Jail is the particle radius, and r is radial distance |14]. Using the weight functions, 
the hard sphere Mayer function can be written as a sum of convolution integrals, 

-fijir) = waiRi) * w-i{Rj) + wi[Ri) * wziRj) + wziRi) * wiiRj) + wsiRO * woiRj), (3.3) 
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where the asterisk denotes the convolution of two functions hi[r) and hzir), which is defined as {hi * 
h2){\r-r'\) - f dx/ji(r-x)^2(r' -x). The spatial arguments of the weight functions have been omitted 
in the notation of the right hand side of <3.3> . The direct space expressions of the weight functions are 
wo{R,r)^-S"{R-r)/{8n) + S'{R-r)/{27ir), wi{R,r) ^ S'{R- r)/[87i), W2[R,r) ^ 5{R- r), and W3[R,r) = 
Q{R- r), where 5 (•) is the Dirac distribution, and ©(•) is the Heaviside step function. 

The appearance of particular combinations of products in l l3.3t can be based on dimensional anal- 
ysis. In order to see this, note that the weight functions Wv are objects with dimensions of (length)^"^, 
hence wo,wi,W2, and carry dimensions of (length)"^, (length)"^, (length)"^, and (length)", respec- 
tively. Each of the products in 1 13. 3> has the dimension (length)"^, which cancels the (length)^ which is 
due to the convolution integral and hence yields a dimensionless Mayer function. 

This property can be formaUzed 133,1 in order to estabhsh more mathematical structure. One can 
re-write <3.3) as 

-fij{r-r']^ [ ^ M"wv(J?i,r-x)w,(J?y,r'-x), (3.4) 

J V,T=0 

where the coefficients M^'' are chosen in such a way that only the "allowed" combinations of weight 
functions contribute. Hence, for most of the index combinations vt, the coefficients vanish, AP''^ - 0. 
Those that contribute in <3.3> possess the index combinations vt = 03, 12,21,31, hence the non-vanishing 
coefficients are M^^ - M^^ - M^^ - M'^" - 1. This procedure amounts to (only) a formalization of 
The crucial step is to view the coefficients M^'^ as the elements of a metric M which operates in the space 
of weight functions w(i?,r). Here, a vector of weight functions w is defined by its components Wv{R,r), 
hence w - {wq, wi, W2, W3). The length scale R acts as a parameter. As an aside, one should not confuse the 
vector w with Rosenfeld's vectorial weight functions W[,i (r) and W[,2 (r) 1 13], which have a very different 
origin, rooted in the geometry of the sphere in three-dimensional space. The w(i?,r), on the other hand, 
are elements of a four-dimensional vector space; their index runs from to 3. 

The metric M can be represented as a matrix, which is akin to a mirrored 4x4 unit matrix in that its 
only non-vanishing elements are unities on the counter diagonal. 



/' 1 A 
10 
10 
10 



(3.5) 



This formalization allows us to write <3.4) [and hence 1 13. 3H more concisely as follows: 

fij = j dxw[Ri)-M-w'[Rj), (3.6) 

where the superscript t indicates matrix transposition, and the spatial arguments have been omitted in 
the notation; these are the same as in ( I3.4> . 

In order to illustrate the framework, we display corresponding diagrams for binary hard sphere mix- 
tures in figurejl] The deconvolution of the Mayer bond /n (r) between particles of species 1 is displayed 
in figure[l](a). The length of the w bonds indicates the magnitude of the argument R, in this case Ri. The 
kink is located at position x, which is the integration variable in the convolution integral I l3.4t . Hence, in 
this and in the following diagrams, the position of a kink (or junction to be introduced below) is integrated 
over. The open circles indicate fixed positions r and r', i.e., the arguments in /n (|r - r'|). One commonly 
refers to these as root points. Recall that multiplying each root point by the one-body density distribution 
and integrating over its position yields (up to a factor of TI2) the ij -II contribution to the exact low 
density limit of the excess free energy functional <3.2) . Figure[l](b) gives the corresponding diagram for 
species 2, here taken to be of a larger size (hence R2 > Ri). The cross species diagram, ij = 12, is shown 
in figure[T](c). From the diagrammatic representation it is clear that the total length of the diagram, i.e., 
the range over which fnir) is non-zero, is fully determined by the accumulated length of the two arms, 
i?i and J?2- Hence, the cross species interaction distance is CT12 - Ri+ R2, a case to which one refers to as 
an additive hard sphere mixture. 
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Figure 1. Diagrammatic representation of the deconvolution of the Mayer bond fijir) between species 
i,j = 1,2 for different types of hard sphere mbctures. The vector of species-specific weight functions 
Wv{Ri,r], where v = 0,1,2,3, is denoted by w. Here, the species index ; = 1,2 is given at the root point 
(open circle). Tensor contraction over pairs of the geometric (Greek) index with the metric M is implied, 
where M is a "reversed unit matrix", i.e., M^^ = Mvt = lifv + T = 3 and otherwise (3.5) . Root points 
(open circles) are fixed spatial positions. Examples are shown for the KR deconvolution for species 11 (a), 
22 (b), and 12 (c) in an additive mixture. For non-additive mixtures, the IDce-species expressions (a) and 
(b) still hold, but the cross-species expressions also possess the kernel matrix K (thick line), as defined 
in (3.9) and (3.10) with index-raised elements K^'^ [d, r] as an additional type of bond. This appears as a 
single bond (d) in the binary non-additive functional, and as a double link in the ternary non-additive 
functional (e), here tensor-contracted by the metric M. This allows one to control the range of fuir) 
arbitrarily, i.e., independently of Ri and i?2. 



Additivity is a fundamental property of the space of weight functions, and is indeed more general than 
outlined so far. Consider a given vector of weight functions w(i?,r), where the length scale R is fixed. We 
aim at finding a linear operation K{d), where d is a length scale, that jdelds 



w^(i?') = K(rf) *w*(i?). 



(3.7) 



where * represents the combined operation of spatial convolution (*) and dot product (■) in the four- 
vector space. Hence, in a less concise but more explicit way I |3.7> can be written componentwise as 



WviR',r-r'] = f dx' Kf-d,r-x]Wr{R,r' -X), 

J v=0 



(3.8) 



where the components of K are denoted by Kl and the two indices v and t run from to 3. That such a K 
exists is a non-trivial matter, and forms the heart of the binary non-additive hard sphere functional fi^, 
see the detailed investigation of the properties of K in 1 33]. The existence and properties are most eas- 
ily demonstrated in Fourier space representation, where the convolutions become mere products. The 
Fourier transform of the "kernel matrix" K possesses a representation as the matrix exponential 



K(i?, q) = exp {RG) , 



where the "generator" is 



G: 



I 
1 








&n 





1 









As an almost trivial consequence of 1 13. 9> . the relationship 

K(i? + R', q) = K{R, q) ■ K(i?', q] 



(3.9) 



(3.10) 



(3.11) 
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holds, simply due to the exponential satisfying exp[{R + R']G) - exp(i?G)exp(i?'G). Correspondingly, in 
real space 

K[R + R',r)^K{R) *K[R'). (3.12) 

Two features of this expression are welcome. One is that the the Kierlik-Rosinberg weight functions 
Wv{R, r) constitute four of the sixteen (4 x 4) components of K{r,R). Hence, the Kierlik-Rosinberg weight 
functions are "automatically" generated when starting with | |3.9> and 1 13. 10> . The second feature is that 
all further components K^iR, r) share their degree of non-locality with the weight functions, i.e., all com- 
ponents of K vanish for distances r > R. This becomes a crucial property when one uses these objects 
as more general convolution kernels in the construction of third and higher orders (in density) in the 
excess free energy functional. Before we describe such developments we return to the problem of a gen- 
eral binary hard sphere mixtures, in which the cross interaction is non-additive, i.e., au ^ Ri+ R2 = 
(an + a22)l2. 

Equation <3.12t allows us to use two convolutions rather than only one in order to represent the cross 
Mayer function for a non-additive binary mixture as 

-fuir] = w(i?i) * K(d) * w*(i?2). (3.13) 

Here, the length scale d is due to the non-additivity, and satisfies a 12 - Ri + d + R2. Figure[T](d) shows the 
corresponding diagrams, where K{d) constitutes an additional type of bond. The length of the bond can 
be adjusted, via changing d, in order to generate the (given) range of fuir). In the diagram, the position 
of each kink is integrated over. These integrals correspond to the convolution integrals in <3.13> . The low 
density expansion of the binary non-additive hard sphere functional |30] contains the same interspecies 
diagrams as the additive Kierlik-Rosinberg functional, as shown in figure [l] (a) and (b), and contains the 
diagram shown in figure[T](d) for the cross species interaction. 

Using the property 13.12> . one can go beyond <3.13> and represent the single K matrix via further 
deconvolution as a (convolution) product of two K matrices, with appropriately chosen length scales di 
and d2, so that ai2 - Ri + di + d2 + R2, and hence 

-/i2(r) =w(J?i) * K(di) * K(d2]*w'(R2). (3.14) 

The corresponding diagram is shown in figure[T](e). This constitutes the cross species low-density limit 
of the ternary non-additive hard sphere functional fsi*]. For ternary mixtures, the three cross species 
lengthscales cri2,0'i3,cr23 are decomposed as 

aij - Ri + di + dj + Rj , (3.15) 

with appropriate values of di,d2 and d-^, which can be uniquely determined for a mixture with three 
components [31]. 

We turn to the third-virial level. The exact contribution to the excess free energy functional is 

L [dr fdr' [ dv"pi[v)pj{v')pk{r")fij{v-r')fj^{v'-v")fik{r-v"). (3.16) 

Due to the pairwise coupling of the three integration variables, one commonly refers to 1 13. 16> as (the 
sum of) triangle diagrams. FMT fails to generate this exact expression , but yields a very reasonable 
approximation to it. Apart from going from Mayer bonds as convolution kernels to weight function bonds, 
the crucial step is a "re-wiring" or topological change in the structure of the diagrams. Avoiding 
the loop in the triangle diagram (a problem that gets more severe with an increasing order in density 
in the virial expansion, see below), the FMT diagrams connect the (three) weight function bonds to a 
common central space integral. A tree-like (or star-like) topology results, see figure [2] Writing out the 
space integrals explicitly yields 

-/,;(r-r')/jfc(r'-r")/;-fc(r-r")== ^ [dxWv(iJi,r-x) w,(i?j,r' -x) u/^CJ^fc-r" -x). (3.17) 

V,T,K = 
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From an algebraic point of view, the operation on the right hand side of i3.17) is similar to a triple scalar 
product of w(i?j), w(i?j), and w(i?j;). The third-rank tensor J that generates this operation can be de- 
rived (gjl from the sole requirement that the result, after the x-integral in 1 13. 17> . does not contain un- 
physical divergences (more precisely that the corresponding order in the partial bulk direct correlation 
function is finite everywhere). The tensor J is symmetric in aU its indices; its non-vanishing components 
are = /^^^ = l, and = I /{An). 




Figure 2. Diagrams of the FMT for additive liard sphere mixtures at the third virial level. On the right 
hand side, the exact third order diagrams are shown for a mixture of two species, 1 and 2, where the 
bonds are Mayer functions fij {r) between species i,j = 1,2. The left hand sides show the corresponding 
FMT approximation. Weight function vectors w for species 1 and 2 are connected to a central junction. 
The position of the central junction is integrated over, and the geometric index of the weight functions is 
tensor-contracted with the third-rank "junction tensor" J. 




Figure 3. Same as figure[2] but for binary non-additive hard sphere mixtures. The diagrams that involve 
either species 1 (a) or species 2 (d) exclusively are the same as those for the additive mixture, compare to 
figure|2](a) and (d). The mixed FMT diagrams (b) and (c) feature the kernel matrix K as an additional type 
of bond. K is a second rank tensor, which is contracted via an M metric to a w bond and, on its other side 
to the third-rank junction tensor J. This aUows one to control the range of the cross Mayer bond /i2('') 
freely. The w bonds are specific for each species (indicated by the species index 1,2 at each root point). 
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Figure[2]illustrates these diagrams for the case of an additive binary mixture. Note that, as on the sec- 
ond virial level, all bond lengths are pre-determined, once the particle radii i?i and R2 have been chosen. 
These are uniquely determined from 7?, = an 12, leaving no room to go beyond the additive case. See (2^ 
for an in-depth discussion of the so-called "lost cases" of FMT that come from the intrinsic differences 
between the triangle and the three-arm star diagrams. In short, when stretching all bonds maximally, the 
exact diagram is larger than the FMT approximation; see |20]. Nevertheless, the different lengths of the 
FMT arms approximate the equilateral cases ijk =111, figure [2] (a) and ijk - 222, figure [2] (d), as well 
as simultaneously the isosceles cases ijk - 112, figure[2](b) and ijk - 122, figure[2](c). For non-additive 
mixtures, however, this is is insufficient for approximating all triangle diagrams. 

For non-additive mixtures, K can be used in order to replace one of the weight functions w(i?) in <3.17) 
by K{d) * w{R). The resulting diagrams are shown in figure[3l The intra-species contributions, ijk =111 
in figure [3] (a) and ijk - 222 in figure [3] (d), are unchanged as compared to the additive case shown in 
figure|2] The inter-species diagrams differ by an additional K bond that acts as a spacer between the arm 
of the "minority"-species (i.e., the one that appears only once, not twice) and the central space integral. 
Note that both "majority" arms directly connect to the center. This applies to ijk - 112, as shown in 
figure [3] (b), and to ijk -\22 as shown in figure [3] (c). 

Although the FMT expressions for third-order diagrams are approximations, they possess one impor- 
tant property which is most clearly analysed when considering the low-density limit of the (partial) bulk 
fluid two-body direct correlation functions, Cjj(r), obtained as second functional derivatives. 



Ci/(|r-r P - 



(3.18) 

p;=const 



Observing (Sj) and | [3A6) 



Cij{r) = /oC-) +EPfcc;,-fc(r) + 0(p?), (3.19) 

k 



where the contribution linear in densities is obtained from the triangle diagram by multiplying one of 
the root points with the bulk density and integrating its position over space, i.e., turning the root point 
of species k into a density field point. The field is constant in the homogeneous bulk. A selection of the 
corresponding integrals is displayed in figure |4] The FMT results give the exact result for c!^.j.(r) for 
all combinations of species. This applies to aU versions of the theory, whether additive (isl binary 
non-additive ISOll . or ternary non-additive (illl . and is a first indication of successful description of bulk 
structure. 

The construction of the additive hard sphere FMT free energy functional can be based on a diagram- 
matic series using the assumption that all diagrams of higher than third order in density are also of 
star-like shape 1241. The weight functions w that constitute the arms are connected to a single central 
space integral. In order to connect p arms of the p-th order in density, a tensor of rank p is required. 
One can show |24] that a recursive relation holds and that these object can be obtained by tensor con- 
traction of an appropriate number {p-2) of third-rank tensors J. As an illustration, at fourth order 
jvv TT _ /'^TT ^ ^j^gj-e t];^e first tensor on the right hand side has one index lowered via contrac- 

tion with the metric, J^^' - Z^/^g r^'^' ■ Together with the coefficients ll[p[p - 1)], which are taken 
from the zero-dimensional properties of the system |24], the resulting series is displayed in figure[5] The 
series can be explicitly summed, and yields, in three spatial dimensions, the Kierlik-Rosinberg form of 
FMT. Each arm together with its filled field point constitutes a "weighted density" (I3lll4il . 



Wr(x) = J drp(r)WT(r-x). (3.20) 

The diagrams that constitute the FMTs for binary and ternary non-additive hard sphere mixture, 
possess different topology; a summary is displayed in figure [6] The one-component hard sphere FMT 
functional, shown in figure [6] (a) consists of stars with arms of the same length R. In one dimension 
(hard rods on a line), this is equivalent to Percus' exact functional 1 26]. The binary version possesses two 
different types of arms corresponding to the small (here species 1) and large (species 2) component, see 
figure [6] (b). The one-dimensional version of the functional for non-additive hard core mixtures fs^l is 
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Figure 4. When one field point is integrated over (with no or a constant density) the approximate FMT 
diagrams give the exact result. Here, a selection of the relevant diagrams is shown. The open symbols 
are root points, i.e., fixed position arguments. The filled symbols are field points that are integrated over. 
Representative examples are shown: the diagram in (a) is for species ijk =111, (b) and (c) are for species 
112, where in (b) a 1-field point and in (c) a 2-field point is integrated over. 




Figure 5. Diagrammatic series of the Kierlik-Rosinberg form of the additive fundamental measure free 
energy functional. The field points (filled circles) represent one-body densities and are integrated over. 
The bonds are fundamental measure weight functions that are joined by an M metric (second order in 
density), or by junction tensor J of fc-th rank (fc-th order in density). The position of the inner junction 
(center of the star) is integrated over. This corresponds to the "outer" integral / dx over the free energy 
density (^{rir (x)) in the standard, weighted density formulation with weighted densities tir (x). The scalar 
coefficient of the fc-th order diagram is 1)]. 

an approximation that differs from the exact solution fss']. The binary non-additive FMT has two central 
junctions, where the arms of each of them belong to one species only (the small species 1 belongs to the 
left center, and the large species 2 to the right center). The two centers are connected via a K bond of 
appropriate length d, so that (T12 = Ri + d + R2, see figure [6] (c). Finally, figure [6] (d) displays a typical 
diagram of the ternary non-additive FMT. Here, there are three subcenters, one for each species, with 
corresponding arms. The subcenters are connected via K bonds to one global center. The positions of all 
junctions are integrated over. A full account of the ternary FMT functional can be found in 131,1 . 

3.2. Non-additive hard sphere fluids in bulk and at interfaces 

We give a brief summary of applications of the FMT for binary non-additive hard sphere mixtures. 
Based on the FMT for this system |3oj,[3l|| but also on Monte Carlo computer simulations, a range of 
physical phenomena were investigated. The bulk fluid structure on the two-body level is described very 
satisfactorily by the theory, as compared to simulation data. This applies to the results for the partial 
pair correlation functions, gij(r), when obtained through the Ornstein-Zernike route using the partial 
direct correlations functions, c, ^ (r), obtained as the second functional derivative d3.18> of the free energy 
functional. This is a quite severe test of the theory, because both the derivative and the Ornstein-Zernike 
equations constitute involved operations. The general performance of the theory is very satisfactory. 
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Figure 6. Examples of the topology of the higher (than third) order diagrams that constitute the respective 
FMT free energy functional. Shown are topologies for the pure system (a), additive binary mixtures (b), 
non-additive binary mixtures (c), and non-additive ternary mixtures (d). Contributions of sixth order in 
density are shown in (a)-(c), and of ninth order in density in (d). 

although the core condition guir < oij) = is only approximately satisfied. This can be rectified with 
Percus' test-particle method [36], where the grand potential is minimized in the presence of an external 
potential that is equal to the inter-particle interaction. Test-particle results for the gijir) reproduce the 
simulation data very well | 37], as long as the system is away from the "decoupling case" of negative non- 
additivity, so that (T12 « (crn -1- (T22)/2. It is important to note that the theory does not yield unphysical 
artifacts in test-particle results for partial pair correlation functions (sTlli^ 

The theory yields analytic expressions for the partial direct correlation functions (both in real space 
and in Fourier space), and hence, via the Ornstein-Zernike relation also analytic expressions for the par- 
tial structure factors Sijiq). This makes it very convenient to carry out an analysis of the asymptotic, 
large distance decay of bulk pair correlation functions from pole analysis of the Sijiq) in the plane of 
complex wavevectors q |39, 40], and to relate these results to the decay of one-body density profiles in 
inhomogeneous situations, such as at interfaces filll . 

For positive non-additivity (T12 > (crn + CT22)/2, which is sufficiently large, the system displays fluid- 
fluid phase separation into two fluid phases with different chemical compositions. The theory gives good 
account of the location of the fluid-fluid demixing binodal in the plane of partial packing fractions of 
the two species (s^lill]. A wealth of interesting interfacial phenomena results as a consequence of the 
bulk fluid demixing, and we refer the reader to the original papers on fluid demixing, asymptotic decay 
of correlations and free fluid interfaces |41], first-order layering and critical wetting transitions in non- 
additive hard sphere mixtures [42] , and capillary condensation of non additive hard sphere mixtures in 
planar confinement l43ll . 

4. Conclusions 

We have described a range of recent developments and applications of classical density functional 
theory in the study of bulk and interfacial properties of liquids. Starting with a reassessment of the un- 
derlying variational principle, we have laid out how to use Levy's constrained search method in order 
to define the free-energy functional. Levy's method provides us with an explicit expression for the free- 
energy functional \2.2\ . We showed that the concept can be generalized in order to define an internal- 
energy functional, which possesses the one-body density distribution and a local entropy distribution as 
trial fields. A dynamical theory built on the internal-energy functional can be found in |12]. Using Levy's 
method, the definition of the functional is explicitly independent of the external potential, which consti- 
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tutes a conceptual advantage. However, the minimization in the function space of many-body probability 
distributions cannot be in practice carried out for a realistic model Hamiltonian. Hence, one has to rely 
on approximations for the functional, as is common in DFT. We refer the reader to Jl2|] for a description of 
various approximate internal-energy functionals, which have been obtained from Legendre transform- 
ing the corresponding approximation for the Helmholtz free-energy functional. 

We have laid out the basic ideas underlying the recent progress in formalizing the mathematical struc- 
ture of FMT. This includes the tensorial nature, in that the "geometric" indices of the Kierlik-Rosinberg 
weight functions play the role of tensorial indices. A diagrammatic notation helps to clarify the non-local 
nature of the excess free energy functional. We have given an overview of applications of the binary 
non-additive hard sphere functional to a range of phenomena. 

The above work extends the previous efforts that were primarily based on the intimate connection of 
the properties of the free energy functional under dimensional crossover, i.e., the result of applying the 
functional to density distributions that correspond to an extreme confinement in one or more spatial di- 
rections, made it possible to generalize FMT to the Asakura-Oosawa-Vrij model [44, 45] of coUoid-polymer 
mixtures (4^, the Widom-RowUnson model 0], and penetrable spheres that interact with a constant 
repulsive plateau fS]. These models have in common that their zero-dimensional properties, i.e., the 
statistical mechanics of a cavity of the size of a single particle, carries still some of the essentials of the 
true three-dimensional problem. Hence, the zero-dimensional problem, which can be solved exactly (or 
with good approximations (aS]) in the above cases, is sufficient as a central modification over the FMT 
for hard spheres, in order to obtain a reliable free energy functional. See e.g., 149-331 for applications to 
confined model colloid-polymer mixtures. 

In future work it would be interesting to apply the dynamical test particle limit fsJ.'ss'] to non-additive 
hard sphere mixtures in order to gain a better understanding of the dynamical behaviour of such mix- 
tures. A further important problem is to re-consider formulating an FMT for soft sphere models ISajsgll 
in the light of the diagrammatic and tensorial structure. Work along these lines is in progress [60]. It 
would also be interesting to investigate the implications of Levy's method for the statistical mechanics of 
quenched-annealed fluid mixtures, where DFT was obtained both via the replica trick i6lil64il . and via a 
first-principles derivation following the Mermin-Evans arguments Issll . 
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HoBi AOC5irHeHH5i b icnaci/mHiii TeopiY ())yHKL4ioHa;iy rycTi/iHi/i: 
(|>yHKL4ioHa/i BHyTpiiuHboY eHepriT i A>cirpaMHa CTpyKTypa TeopiV 
(|>yHAaMeHTa/ibHoY Mipi/i 

M. LUMifliGEl, M. Bypricp], B.C.B. AsaHAap^, T. /IflMia/ibP, n. TonKiHcP 

Biflfli/ieHH^i TeopeTWHHoT c|)i3HKM II, iHCTHTyT (fiisnKM, YHiBepcuTeT BaiJipoi^Ty, 95440 BaiJipoiiiT, HiMenHUHa 
tt)i3M4Ha ;ia5opaTopi5i fJ. Bi/i/iic, VHiBepcwTeT BpucTO/iK), BpncTO/ib BS8 1TL, Be/inKa Bpi/iraHi^i 
Katjjeflpa (fiisnKi/i, flepxasHMM yHiBepci/ireT flxoK'^Kapii/i, flxoK'^Kapra, iHflOHesifl 

3po6;ieHO or/iflfl fleKi/ibKOx HeflasHix npaL4b 3 leopiT (fiyHKi^iOHa/iy rycTHHH fl/ia K;iacM4HMX HeoflHopiflHux pi- 
fli/iH. Mn noKasycMO nmM huhom Merofl /leei o6Me>KeHoro noiuyKy MO>Ke 6yTH BUKopucraHnPi fl/ia BUBefleHHfl 
BapiaL(ii^Horo npMHL4Mny, akmm /lexnTb b ochobI reopii c|)yHKL4iOHa;iy rycTHHH. HepeBaroio L4boro Merofly e re, 
mo Bi/ibHa eHeprifl re/ibMro/ibma c|)yHKL4iOHa;i npo6HoT OflHonacTUHKOBOi rycTHHH saflacTbca y ^iBHOMy bh- 
Mflfli 6e3 BiflHOCHO flo BOBHiiuHboro n0TeHL4ia;iy, i^e e y BanaflKy CTaHflapiHoro flOBefleHHfl MepiviiHa-EBaHca 
HepeB reductio ad absurdum. Mm noKasycMO AKysara/ibHUTW niflxiflfl/ia Toro, 1406 BnpasHTW BHyipiLUHK) enep- 
rlK) y BMMflfli c(>yHKL4iOHa;iiB posnofli/iiB OflHonacTUHKOBoT rycTHHH i /lOKa/ibHoT eHiponiT. Tyr /lOKa/ibHuPi xImI- 
4Hi/ii^ noTeHL^ia/i i leMneparypa b o6'eMi Biflirparoib po/ib MHO>KHMKiB /larpanxa b pIbhahhI EM/iepa-ZlarpaHxa 
fl/ifl MiHiMi3aL4M ((jyHKL^iOHayiy. Ha5;iwxeHHfl flyia c|)yHKL4iOHa;iy Bi/ibHoT eHepriT, noKasaHO fliarpaMHy crpyKiy- 
py (fiyHKL^iOHa/iy rycTWHW c|)yHflaMeHTa;ibHOi Mipn PoseHcfieyibfla fl/ia cyMiiui TBepflwx cc|)ep. OnwcaHO HeflaBHi 
ysara/ibHeHHfl fl/ia 6iHapHwx i noTpiwHux cyMiujePi TBepfli/ix C(()ep, akI rpyHiyioTbCfl Ha CKa/iapnux SBaxyBa/ib- 
Hi/ix cfiyHKL^mx Kiep;iiKa-Po3iH6epra. 

KjiiOMOBi cyiOBa: Teopia (pyHKi^ionany rycmHH, reopeMa F oreH6epra-KoHa, (pyHKi^ionan PosBHCpejibfla 
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